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The Debye- Waller factor has been calculated for stabilized <5-phase plutonium with 5% aluminum. A quasi- 
harmonic Born-von Karman force model with temperature dependent phonon frequencies was used to calculate 
the mean-square thermal atomic displacement from absolute zero to 800 K. Implementation of the observed 
anomalous softening of the long wavelength phonons with increasing temperature cannot account for the soft- 
ening of the measured thermal parameter at high temperatures nor for its rather high value at low temperatures. 
The implications for diffraction measurements on 5-phase stabilized plutonium alloys are discussed. 



INTRODUCTION 

Since the measurement of the phonon dispersion curves along the high symmetry axes of small crystallites of (5-phase 
^■^^PUo.98Gao.2 ma, as well as measurements of the phonon density of states and elastic moduli on polycrystalline sam- 
ples of ^''^Puo.gsAlo.os fl, it has become possible to construct phenomenological phonon models which accurately reproduce 
the measured specific heat all the way up to 300 K 101 • Unfortunately, the vibrational frequencies or density of states were 
determined only at a few temperatures. However, from measurements of the elastic shear and bulk moduli it is known that the 
long wavelength frequencies of (5-stabilized Pu exhibit an anomalously large softening, which is nearly linear from absolute zero 
to the melting temperature (see Fig. 1) JHIHSl- 

Here, I address the question of how accurately such a quasi-harmonic phonon model, with temperature dependent frequencies, 
can reproduce the mean-square thermal atomic displacement parameters measured by Lawson et al. HJj, 3t Qi- The working 
assumption is that all force constants, i.e, all frequencies, can be scaled by the same temperature dependent function as the 
observed elastic moduli or long wavelength phonons. This procedure worked very well for the calculation of the specific heat of 
242puo.95Alo.o5 between 2 K and 300 K lAifll. 



RESULTS AND DISCUSSION 



The exponent of the Debye-Waller factor e^^^'^'^^ (mean-square thermal atomic displacement parameter) is given by the 
average of the square of the scalar product of the relative atomic displacement u from its average position with the neutron 
scattering vector q, for example, see Ref. m for details: 2H/(q) = ((q-u)^). In the case of a monatomic cubic crystal this 
simplifies to the standard result for the thermal parameter 2W{q) = 2Woq^, where 2Wo — {u^) = (Uy) = (ul). The isotropic 



thermal parameter is given by 111 U 112 1 
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Here M is the atomic mass, uj is the phonon frequency and the phonon density of states is normalized to satisfy dcjN{uj) — 1. 
At temperatures much larger than the Debye temperature 8 d the Debye- Waller temperature Qdw can easily be extracted from 
the asymptotic form of 2Wq{T) oc T/8|)^. For a harmonic crystal, Qdw is a constant, while Lawson et al. have found that 
for most light actinides 8 dw acquires a temperature dependence. 



Figure 1 shows the temperature dependence of the shear modulus G of (5-phase stabilized plutonium alloys. Although there is 
variation in the magnitude among those different alloys, their slopes are surprisingly similar. Using the temperature dependence 
of G reported for ^^^Puo.gsAlo.os and the 3NN (third next-neighbor) Born-von Karman force model discussed in Ref. with 
parameters given in Table H] I calculated the mean-square thermal atomic displacement parameter 2Wq{T) as defined in Eq. 1. 
The phonon dispersions were scaled to give the following long wavelength elastic coefficients Cij{T) — /(T)Cij (300K), with 
scaling function /(T) = 1.26 - 8.67 • IQ-^T/K, and Cii(300K) = 34.0 GPa, Ci2(300K) = 24.9GPa, and C44(300K) = 31.4 
GPa. Note that these single crystal moduli are consistent with the measured polycrystalline shear, G, and bulk, B, moduli. In 
Figure 2, 1 compare the theoretical 2 Wq (T), with no fit parameters, with the measured mean-square thermal atomic displacement 
parameters by Lawson et al. 101 ■ It is worth to point out that Lawson et al. measured the low temperature data in a cryostat 
on HIPD at LANSCE (Lujan Neutron Scattering Center at Los Alamos National Laboratory) and the high temperature data 
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FIG. 1: Temperature dependence of the shear modulus G of various polycrystalline 5-phase Pu samples (1 wt-% Ga ~ 3.35 at-% Ga) fllBS]- 
The dashed lines are linear fits to the data of Refs. 13li5i]. Note a very similar temperature dependence is reported for the bulk modulus B. 



TABLE I: Bom-von Karman force constants ^f^- for an effective monatomic /cc crystal Puo.95Alo.05 at 300 K with lattice parameter 
ao = 4.59A. For simplicity no correction was made to account for the presence of 5% Al. 
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8.87885 


-1.20410 


11.92267 
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-1.29065 


0.52554 
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-0.93343 


0.03871 


-0.48034 


-0.78982 



in a furnace on NPD at LANSCE. For comparison 1 show the thermal parameters with T-independent force constants, i.e., T- 
independent frequencies, which were arbitrarily determined at T = 400 K for better agreement with experiment. The calculation 
with the frozen-in force constants is in overall better agreement with experiment than the one with T-dependent force constants. 
This agreement is most likely fortuitous, but emphasizes the weak deviation of the experimental 2T4^o from linearity. 

It is surprising that the realistic T-dependent phonon model shows the largest deviations from experiment at low and high 
temperatures. At high temperatures the Debye-Waller phonon moment, Qdw, which coiTesponds to 2Wq{T), samples mostly 
the long wavelength frequencies l/oj'^) and therefore agreement should improve. Both model calculations show significant 
deviations from the experimental 21^0 (T) at low and high temperatures, while agreement is good between 180 K and 500 K. 
Lawson and coworkers argued that the thermal parameter can only be determined up to an additive constant by assigning a large 
offset of 40-80 pm^ to instrumental and systematic uncertainties in the Rietveld refinement of their neutron powder diffraction 
patterns. In later work Lawson and coworkers [8] avoided discussing the nature of this large offset all together. This offset in 
the thermal parameter is much larger than typically seen in diffraction measurements like fee Ni, fee Ce and GaAs ifisi [l4ll 
when comparing experiment with simple phonon models. Although 7-Ce and S-Pn show many similarities among fee crystals 
the large measured low-T value for 2Wo{T) cannot be explained simply by scaling the hypothetical absolute zero temperature 
thermal parameter of 7-Ce (2Wo(0) ~ 25 pm^) with the ratio of the atomic masses of Ce over Pu (~ 140/242) [13^, ,14fl . 
Here, we obtain for 6-Pu the theoretical absolute zero temperature value 2Wo(0) ~ 13 pm^, which is much smaller than the 
experimental value of 2Vl/o(10K) « 49 pm^. Therefore, the large offset observed by Lawson might be due to a combination 
of systematic and instrumental errors as well as static disorder induced broadening in the static structure function, which can 
lead to an overestimate of the mean-square thermal atomic displacement parameter in a Rietveld refinement. This interpretation 
is consistent with the more recent observation by Lawson and coworkers [15] that there is significant contribution of diffuse 
scattering in neutron powder diffraction patterns of Puo.98Gao.02- Note that the low-temperature and long-wavelength phonon 
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FIG. 2: The mean-square thermal atomic displacement parameter 2Wo{T) vs. temperature for Pu with 5-at.% Al. A comparison is shown 
between Lawsons's data (circles) [Tj] and (i) a T-dependent (solid line) and (ii) T-independent (dashed line) 3NN Born-von Karman force 
model. For case (ii) the force constants where frozen-in by using the values of model (i) at T = 400 K. Consequently Qdw is constant for 
T > ec ^ 116 Kd. 



model successfully reproduced the low-temperature specific heat measurements fl, thus it is not clear why it should fail to 
explain the low-temperature thermal parameter. On the other side, the deviations seen at high temperatures could be due to 
anharmonicities and the breakdown of the simple scaling function used for modeling the T-dependent frequencies. However, 
the breakdown of the scaling function for the phonon dispersion is unexpected and requires further diffraction studies because it 
means that above ^ 500 K the phonons show no longer the anomalous softening as seen in the elastic moduli (see Fig. 1). 

CONCLUSIONS 

A 3NN Born-von Karman model with force constants frozen in at 400 K fortuitously reproduces the measured Debye-Waller 
factor over a wide temperature range, while a model with more realistic temperature-dependent phonon frequencies overesti- 
mates the mean-square thermal atomic displacement parameter above ~ 500 K. Both models fail to explain the large experi- 
mental low-temperature thermal parameter below 100 K and its large asymptotic zero-temperature value (~ 49 pm^) due to 
zero-point motion. The experimental low-T mean-square thermal atomic displacement parameter is nearly four times bigger 
than the theoretical value or the corresponding value of 7-Ce scaled by the atomic mass ratio. These discrepancies are surpris- 
ing since the same temperature-dependent phonon model was very successful in describing the specific heat of ^^^Puo.gsAlo.os 
between 2 K and 300 K. Further diffraction studies are needed to investigate the origin of the large deviations between theory 
and experiment at low and high temperatures, as well as the large unphysical offsets used in previous neutron powder diffraction 
measurements. 
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